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We investigate the entanglement spectra arising from sharp real-space partitions of the system 
for quantum Hall states. These partitions differ from the previously utilized orbital and particle 
partitions and reveal complementary aspects of the physics of these topologically ordered systems. 
We show, by constructing one to one maps to the particle partition entanglement spectra, that the 
counting of the real-space entanglement spectra levels for different particle number sectors versus 
their angular momentum along the spatial partition boundary is equal to the counting of states 
for the system with a number of (unpinned) bulk quasiholes excitations corresponding to the same 
particle and flux numbers. This proves that, for an ideal model state described by a conformal 
field theory, the real-space entanglement spectra level counting is bounded by the counting of the 
conformal field theory edge modes. This bound is known to be saturated in the thermodynamic 
limit (and at finite sizes for certain states). Numerically analyzing several ideal model states, we 
find that the real-space entanglement spectra indeed display the edge modes dispersion relations 
expected from their corresponding conformal field theories. We also numerically find that the 
real-space entanglement spectra of Coulomb interaction ground states exhibit a series of branches, 
which we relate to the model state and (above an entanglement gap) to its quasiparticle-quasihole 
excitations. We also numerically compute the entanglement entropy for the = 1 integer quantum 
Hall state with real-space partitions and compare against the analytic prediction. We find that the 
entanglement entropy indeed scales linearly with the boundary length for large enough systems, but 
that the attainable system sizes are still too small to provide a reliable extraction of the sub-leading 
topological entanglement entropy term. 

PACS numbers: 03.67.Mn, 05.30.Pr, 73.43.-f 



I. INTRODUCTION 

Quantum information theory provides powerful tools 
and perspectives for understanding and characterizing 
quantum many-body systems. This is particularly true 
for topologically ordered phases of matter. Such systems 
possess profoundly entangled states and cannot be char- 
acterized by local order parameters, so it is natural to an- 
alyze and characterize them through their entanglement 
properties. One of the most basic (quantum information 
theoretic) tools for studying entanglement is entropy. If 
a system in a pure state p is partitioned into two subsys- 
tems A and then their entanglement entropy, i.e. the 
von Neumann entropy 

Sa = Sb = -Tr [pA log Pa] = -Tr [ps log pb] (1) 

of either reduced density matrix pA = Tr^ [p] or pB = 
TiA [p], is the unique measure of entanglement between 
the two subsystems. 

For a real-space partition of a gapped, two-dimensional 
(2D) system, the leading order contribution to the en- 
tanglement entropy obeys a "perimeter law," i.e., it is 
proportional to the length C of the boundary, with a 
non-universal coefficient a, as £ ^ oc^ ^. Kitaev and 
Preskill^ and Levin and Wen^ showed that the sub- 
leading, constant term of the entanglement entropy for 
such a system in its ground state is a universal quan- 
tity, which they called the "topological entanglement en- 
tropy." In particular, the entanglement entropy (for the 



ground state) takes the form Sa = <^C — 77,7, up to terms 
that vanish in the limit £ ^ 00, where n is the number 
of connected components of the boundary of A. Here, 
7 = logP, where P > 1 is a quantity known as the "to- 
tal quantum dimension," which is characteristic of the 
system's topological order and equal to 1 only for trivial 
topological order (i.e., a gapped system with no anyonic 
excitations)^^. 

Quantum Hall systems are, so far, the most studied 
and only experimentally realized topologically ordered 
phases. Most observed Hall plateaus are expected to be 
Abelian quantum Hall states^ with quasiparticles that 
have fractional charge and exchange (braiding) statis- 
tic4^. However, some of the observed plateaus in the 
second LandaiP^HI^, most notably: u = 5/2 and 12/5, 
are expected to host non- Abelian statea^^^^, which pos- 
sess quasiparticles with non- Abelian braiding statisticP^. 
These non- Abelian states have recently received much in- 
terest due to their potential application for topologically- 
protected fault-tolerant quantum computation^*^. 

The entanglement entropy depends on the way the sys- 
tem has been partitioned. Until now, two kinds of bipar- 
tite partitions have been applied to fractional quantum 
Hall (FQH) states: the orbital (or momentum) parti- 
tion^^, for which A consists of the first I a orbitals (on 
the sphere or torus) while B consists of the remaining 
ones and the particle partition^^, for which A consists of 
the first Na particles and B consists of the TV — Na re- 
maining ones. A third partition, the real-space partition. 
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for which A is a spatial region and B is the complement of 
A, has been applied to the integer quantum Hall (IQH) 
state^^, which can be treated analytically. These dif- 
ferent partitions provide access to different properties of 
the state. Using an orbital partition as an approximation 
of a real-space partition (which was argued to be a rea- 
sonable approximation, since the orbitals are Gaussian 
localized), Haque et al?^ have numerically extracted the 
expected value 7 = ^ In 3 for the v = !/?> Laughlin state. 
However, the extraction of the topological entanglement 
entropy (using the orbital partition) in numerical studies 
has proven to be far less reliable for anything but the sim- 
plest model wavefunctions, in particular for the Coulomb 
Hamiltonian ground stated. 

Even though the topological entanglement entropy can 
provide significant information regarding the topological 
order of a system, it is not able to uniquely identify it. In 
Ref.T, Kitaev and Preskill provided an additional heuris- 
tic derivation of the topological entanglement entropy, in 
which they suggested that the reduced density matrix pA 
for a disk-like region A can be used to define a Hamilto- 
nian 
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FIG. 1: (a) The OES for the = 1/3 Laughlin state with 
= 8 and A^A = 4. (b) The OES for the z/ = 2/3 particle- 
hole conjugate of the = 1/3 Laughlin state with = 14 and 
Na — 7. The two spectra are mirror images of each other: 
the v — spectrum contains a chiral mode, moving upward 
toward the left, while the = 2/3 spectrum contains an anti- 
chiral mode, moving upward toward the right. The counter- 
propagating boson mode that should appear in the z/ = 2/3 
state's edge mode spectra does not appear in the OES for the 
same reason the v — 1 integer quantum Hall state exhibits 
trivial OES, i.e., the OES probes only the "interacting" part 
of the wavefunction. (We note that the degeneracy between 
states at L\ = 24,23 for v = 1/3 spectrum is a detail due 
to the fact that, for this particular partition, the two states 
enjoy an L symmetry, which is not present for some other 
partitions nor for the bosonic version of this state.) 



Ua = - log PA (2) 

which, for a topologically ordered system, would be natu- 
ral to regard as the Hamiltonian of a (1 + 1) -dimensional 
conformal field theory (GET) corresponding to the edge 
of the system, at least for the purposes of determining the 
universal contribution to the entanglement entropy. In 
an effort to develop a more discriminating tool for iden- 
tifying the topological order of a system (particularly for 
non-Abelian FQH states), Li and Haldane^^ proposed ex- 
amining the (low- lying portion of the) spectrum of the fic- 
titious Hamiltonian 1-La = — logp^, which they dubbed 
the "entanglement spectrum" (ES), versus the quantum 
numbers of the available symmetries in the problem. 

For both the Moore-Read (MR) model stat^and the 
ground state of the u = 5/2 Goulomb interaction, Li and 
Haldane^^ plotted the ES versus the quantum num- 
ber of angular momentum along the z-axis for subsystem 
A. For computational expediency, they used the orbital 
partition, assuming that, since the orbital basis states are 
Gaussian localized, this would give a good approxima- 
tion to the real-space partition. They found that the ES 
matched the expected MR GET [i.e., Ising x U(l)] edge 
mode counting at the lowest pseudo-energies. Given this 
numerical observation, they suggested (along the lines 
of Kitaev and Preskill's statement) that the counting 
of the low-lying ES levels would match the GET edge 
mode counting for all quantum Hall states (up to some 
limit, which grows with the system size). This obser- 
vation is even more remarkable if one realizes that for 
generic states the number of entanglement levels should 
saturate the Hilbert space dimension, a number expo- 
nentially larger than the number of GFT edge modes. 
Renormalization group arguments supporting this con- 
jecture have been given in Ref. 31 , whereas Ref. 32 proved 
that the number of levels of the orbital partition ES is 



bounded from above by the number of the GFT edge 
modes, thereby proving a large part of the conjecture. In 
the absence of accidental reductions of the rank of the 
density matrix, this bound should be saturated. 

However, since the orbitals have non-zero support 
across the entire sphere (or plane), the orbital partition 
entanglement spectra (OES) will not precisely match the 
real-space partition entanglement spectra (RES). A sim- 
ple example of this mismatch is provided by the u = 1 
IQH state, for which the ground-state wavefunction is 
the fully filled n = Landau level. An easy inspection 
reveals that the ground state p is a product state in the 
orbital basis and, hence, pA is a pure state for any or- 
bital bi-partition (where all Na orbitals in subsystem A 
are filled) of this system. Gonsequently, the ES (and, 
hence, the entanglement entropy) is trivial, consisting of 
exactly one level, rather than exhibiting the edge GFT. 
On the other hand, the edge mode counting is expected 
to be that of a U(l) GFT (i.e., a free chiral boson). The 
RES has b een computed analytically for non-interacting 
state^^^^, which is done by filling up the levels of a 
single-particle spectrum (an easily obtainable quantity). 
For such states, the RES was found to agree with the 
edge GFT, with a linearly scaling entanglement entropy. 
Additional examples of the mismatch between the OES 
and RES are provided by particle-hole conjugation^^ of 
quantum Hall states (such as the v = 2/3 state obtained 
from the = 1/3 Laughlin state by particle- hole conju- 
gation) . The OES of a state is exactly the mirror image, 
up to a constant shift for L^, of that of its particle- hole 
conjugate state, see Fig. [ijfor the u = 1/3 and u = 2/3 
Laughlin state example. The RES is expected (from the 
edge GFT) to be given by a combination of the mirror 
modes with an additional mode propagating counter to 
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these. In a sense, the orbital partition only probes the 
"interacting" part of the state and is not sensitive to the 
presence of the (non-interacting) integer- filled portion of 
a quantum Hall state. 

In this paper, we examine the ES of quantum Hall 
states using real-space partitions, both on the sphere and 
on the disk geometry. The paper is organized as follows: 
In Section [Hj we provide the general formalism for calcu- 
lating r eal-sp ace partition reduced density matrices. In 
Section II A , we specialize the formalism t o (bos onic and 
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we show 



fermionic) quantum Hall states. In Section '. 
that the RES is related to the particle partition entan- 
glement spectrum (PES)^^ through a one to one map. 
For quantum Hall states described by a GET, i.e., with 
trial wavefunctions generated by conformal blocks, this 
proves that the RES level counting is bounded by the 
GET edge mode level counting. Moreover, it proves that 
this bound is saturated in the thermodynamic limit for 
such states and at finite sizes for some states. We expect 
that the bound will also be saturated at finite sizes by all 
such states. (This is supported by our numerical results.) 
In Section [Tin we numerically analyze several important 
examples. In particular, we analyze the v = 1 IQH state, 
the = 1/3 Laughlin state and its v = 2/3 particle- 
hole conjugate state, and the v = 1 bosonic MR state, 
which all have RES that exhibit the structures expected 
from the corresponding GETs. We also analyze the RES 
of the ground-states of realistic Goulomb Hamiltonians 
and find that they exhibit the same low-lying levels as 
the corresponding model states believed to describe their 
universality classes. The RES of ground states of these 



Hamiltonians also exhibit several other branches, besides 
that of the model states', in the higher ES levels. In Sec- 
tion |IV| we relate these higher branches to non-Abelian 
quasielectron-quasihole excitations of the system for the 
case of the v 



5/2 Goulomb state. In Section VI 



we an- 
alyze the RES for the Laughlin model state as we change 
the geometry of the manifold from the disk to the annulus 
to show how the edges of the sample infiuence the RES. 
In Section [V| we analyze the RES of the MR state with 
quasiholes. We show that, when a non-Abelian quasihole 
crosses the partition boundary, the RES exhibits a sharp 
change corresponding to the change of topological sec- 
tors of the GET edge mode. Finally, in Section |VII[ we 
compute the entanglement entropy using the spatial par- 
tition for the V = 1 IQH state to verify the validity of the 
entanglement entropy perimeter law and try to extract 
the topological entanglement entropy. While we find the 
entanglement entropy scales linearly with the length of 
the boundary (for long boundary length), but are unable 
to reliably compute the topological entanglement entropy 
even for this simple case. This shows that determining 
the topological entropy directly from numerics without 
using prior knowledge about the state is not currently 
possible. 

II. REAL-SPACE ENTANGLEMENT 
SPECTRUM 



We consider a pure, N particle state with wavefunction 
^ and a bipartite real-space partition into regions A and 
B. The (normalized) density matrix is 



. , , . ^*(ri,...,rAr)^(r^,...,rV) ... 
p(ri, . . . , Tat; r 1, . . . , r at) = . (3) 

/n^^r,|^(ri,...,r^)|^ 

(We allow for unnormalized wavefunctions, but the density matrices and reduced density matrices throughout this 
paper are properly normalized to have unit trace.) The reduced density matrix p^, obtained by tracing out the 
degrees of freedom in region 5, will commute with 7V^, the number operator of particles in region A, so the density 
matrix is block diagonal in sectors of different Na-, which can be considered separately. In particular, 

N 

PA= ^ Piv^PAr^(ri,...,rAr^;r'i,...,rVA), (4) 

Na=0 

where 



f N\ f f ^ 

P^A = [.r / Tl^^^i TT p(ri,...,rAr;ri,...,rAr) (5) 
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is the probabihty of the reduced system having Na electrons in region A with (normahzed) density matrix 

N 

Ib n d'^^j p(ri,...,rAr^,rAr^+i,...,rAr;r'i,...,rVA,rAr^+i,---,rAr) 

/)Ar^(ri,...,rAr^;r i,...,r AT^) = 



Na N 

Ia n d'^^ils n d^^j p(ri,...,riv;ri,...,rAr) 

i=l J^Na + I 



1 /iv— ^ 



\ AT ) TT c^^rj P(ri,...,rAr^,rAr^+i,...,rAr;r'i,...,r'Ar^,rAr^+i,-.-,rAr) (6) 

where the coordinates of the Na remaining particles are restricted to region A (and the Nb = N — Na particles 
confined to region B have been traced out). We note that if p is a product of (identical) single-particle states, then 
PNa = {na)Pa^^^ — Pa)^~^^, where pa is the probability for any one of the particles to be in region A, is a binomial 
distribution (which becomes a Gaussian distribution as A/" ^ oo). In the rest of the paper, we will work with the 
density matrix pN^ of a specific Na subspace, and leave the coordinates implicit when the context is clear. 

A. Quantum Hall States 

A general bosonic or fermionic quantum Hall wavefunction for N particles in the lowest Landau level with total 
z-axis angular momentum L^^^ along the z-axis can be expressed as linear combinations of Fock states of single 
particle orbitals: ^ = cxA4x. Here A = (Ai, . . . , Aat) is a length N partition of I/tot 5 ^•^•5 ll'^ll = Xli — ^tot 
\i > A^+i. The normalized symmetric monomial (for bosonic systems) or Slater determinant (for fermionic systems) 
that has orbital occupation corresponding to the partition A is given by 

1 ^ 

= mnnmmS ^^"\n/A..(r.), (7) 

m ^ 

where Sn is the permutation group of N objects, e(<j) = (±1)^ is equal to 1 for bosonic systems or the signature 
of the permutation a for fermionic systems, nm(A) is the number of occurrences of the integer m in the partition A, 
corresponding to an occupation of the mth orbital (which are all equal to or 1 for the fermionic case), and 



0m (r) 



V27r2'^m! 



z^e 4 1^1^ plane 



4.iy;2L)! [cos(^/2)e'W2]™ [sin(^/2)e-^/2]^*-'" sphere 



(8) 



is the single particle orbital with eigenvalue m on the plane with coordinates z = x -\- iy (where m = 0, 1, . . .), or 
with eigenvalue m — N(j)/2 of the sphere with polar coordinates (^, ip) and A/'^ fiux quanta through the surface of 
the sphere (where m = 0, 1, . . . , A/'^). The functions Mx form a set of orthonormal free many-body states. 

Using Eq. (|6|, we can derive the expression of the reduced density matrix for the real space partition as a function 
of the coefficients cx'. 



AT 



n / ^'r,^^.o)(ri)'/'V,,(r,). (9) 



To proceed further, we need to specify the spatial regions A and B. While a generic integration domain can certainly 
be used, we want to minimize our computational burden and also be able to compare the resulting ES with that of 
previous approaches. For these reasons, we chose to use domains that are rotationally invariant around the z-axis, 
which preserves the symmetry. The projection on the z-axis of the angular momentum for the particles in region 
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A, L^, was key in identifying the model state structure within the ES of reahstic systenPSl. On the plane, we take 
region A to be the disk-shaped domain centered at the origin with radius R (i.e. the region r < R). On the sphere, we 
take A to be the cap centered at the north pole up to azimuthal angle O (i.e., the region ^ < O). For such domains, 
the integrals of Eq. (|9| will lead to a vanishing contribution unless Xa{j) = ^T{j) ^ ~ ^ 1, . . . , A/". This 

ensures that [pA^L\] = and, thus, pA is also block-diagonal in L\ sectors. To write pN^ more explicitly, we use the 
property 

A^A(ri,...,r^) = ^(m,^) i ^^j^''' n [nrn{^l)^^^^^ >1;.(ri, . . . , r^^)M,(r^^+i, . . . , r^) (10) 



where ijl are length Na partitions, u are length Nb partitions, we define (/[x; u) to be the ordered partition with elements 
/ii, . . ..pnai^Ii ' ' . ,i^Ars, we define e{fi,u) = e{a) for a e Sn such that (A^(i), . . . , A^(Ar)) = (/ii, • • . ,/iAr^,i^i, . . ^,J^Nb) 
for A = {/J^]iy) (i.e. <j~^ is the permutation needed to combine fj, and ly into an ordered, length TV partition), and 
the sum is over all such partitions that can be combined to give the partition A (i.e. /x and i/ are subjected to the 
constraint {/j,; v) = A). We also define 

r(5(m + l,^) plane 
^B(m) = y d\\<i>^(v)\^ =\ (11) 
^ [^cos2(f )(^ + l^^c/) - ^ + 1) sphere 

to be the norm-squared of the orbitals restricted to region 5, where Q(a^x) = r(a,x)/r(a) is the regularized incom- 
plete gamma function and /^^(a, 6) = B(x\ a, h)/B{a^ b) is the regularized incomplete beta function. This allows us to 
write 

PN^ = K,.R^',^M*^M^^ (12) 



We note that this indicates that the RES will weight the different orbitals according to the overlaps of orbitals in the 
traced out region B. 

We can also write this in terms of L\ sectors as 

PNa = p{L\\Na)pna,l^^ (14) 



Nb 



where p{LWNa) is the conditional probability of the re- 
duced system A having angular momentum L^, given 
that it has Na particles, andpAr^,!.^ = p{La\^a)pna is 
the probability that the reduced system A has Na parti- 
cles and angular momentum L\. 



We emphasize that the coordinates are restricted to 
the region A and so one might want to write the reduced 
density matrices in terms of basis states that are normal- 



ized on region A. Defining 



Mt = 



0m (r) 



Na 



(16) 
(17) 



A l[^^-^B{p^)] R^.,. (18) 



it is clear that we can simply substitute these in for their 
corresponding quantities in the above equations. Using 
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Mf^ to define orthonormal basis states \ij,a) for subsys- 
tem A, we can ttius write ttie density matrix operator 
as 

pN^= E K:.K',.WA){t^A\. (19) 

In otiier words, tiie density matrix elements in this basis 
are 

[pNA]fj^',fj^ = X] = (20) 

where on the right hand side we are treating R"^ as a 
matrix with row indices /j. and column indices ly. 

It is worth pointing out a distinction between the or- 
bital and real-space partitions that generally appears in 
the probability distribution p^^ of having Na particles 
in subsystem A. As the (non-interacting) IQH states 
are product states, the probabilities form a bino- 
mial distribution. On the other hand, applying an or- 
bital partition to an IQH state gives = ^Na^N"^^^ 
where Nj^^ is the total number of orbitals comprising 
subsystem A. More generally, one expects the interact- 



ing quantum Hall ground states to exhibit a probabil- 
ity distribution pn^ that is qualitatively similar for the 
real-space partition, i.e., a nearly binomial distribution 
peaked around Na ^ /aN, where /a = i^^^^u^- Simi- 
larly, the probability distribution for the orbital partition 
will be sharply peaked around Na ~ /a^, where here 
Ia — ^orb^^^orb 7 and will drop to exactly zero for values 
of Na beyond some range away from JaN^ as dictated 
by the squeezing properties of the state. 



B. Relation between real- space partition and 
particle partition entanglement spectra 

We can now demonstrate that there is a one to one 
map between the RES and the PES of an arbitrary N 
particle state |^). For a particle partition of the system 
into subsystems A and B comprised of N a and Nb par- 
ticles, respectively, the reduced density matrix operator 
for subsystem A is obtained by removing Nb particles, 
irrespective of their position in space. This is given by 



p^(ri, . . . , tat^; r 1, . . . , r ata) = / H ^^^^ ^^^i' • • • ' ^Na^^Na+i^ • • • , tat; r i, . . . , tn^.tna+i, • • • , tat) (21) 



j=NA+l 



where the domains of the integrals and the coordinates of the remaining particles are the entire physical space of 
the system. We notice, apart from normalizations and the difference of domains of the integrations and remaining 
coordinates, this has a form similar to pna of Eq. (|6|, so it should be clear that the two can be related. 
Using the decomposition of Eq. (10), we can write the wavefunction ^ = ^x^><'^^ 



Applying the particle partition, the reduced density matrix operator's elements are given simply by 
where we have defined the matrix P with elements 



I 



with row indices /j. and column indices ly. 



Comparing this to pN^ in Eq. pQ| ), we can define the 
diagonal matrices Q and 5* with elements 



1 



\ PNa 



^ / j=i 



7 



and 



[S] 



1 



Na 

nil 



(26) 



to give us the relation R"^ = SR = SPQ. Within the 
many-body basis, the calculations of the RES can now 
be done numerically in a way similar to the calculations 
of PES. 

Since Q and S are diagonal with no zeros on the diag- 
onal, i?^, i?, and P have the same rank. Consequently, 
the reduced density matrix operators of the particle 
partition and p^^ of the real-space partition in the Na 
sector have the same rank. It is clear that this also holds 
when one restricts to a particular L\ or sector. How- 
ever, it is worth mentioning that, on the sphere, the total 
angular momentum L\ of subsytem ^4 is a good quantum 
number for the PES, but is not a good quantum num- 
ber for the RES and the OES. This gives a one to one 
mapping between the levels of the RES and the PES. For 
states described by a CFT, the counting of the PES levels 
per momentum sector was shown to equal the counting 
of states for the system with unpinned bulk quasiholes 
corresponding to the same number of particles Na and 
orbitals A/'^ + 1^^. It has been proven for such quantum 
Hall states that this counting of states with quasiholes 
is bounded by the level counting of the corresponding 
CFT edge modes^^. Moreover it has been proven that 
this bound is saturated in the thermodynamic limit for 
all such quantum Hall states and also at finite sizes for 
Laughlin stateJ^. Thus, we have proven that the RES 
level counting of such model states similarly satisfies and 
saturates the bound by the CFT edge mode counting. 



III. EXAMPLES 
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FIG. 2: The RES of the zy = 1 slater state for iV = 11 particles 
with Na — 5. The sphere was partition into two hemispheres 
(e = 7r/2). Inset: The OES for the same state {Na = 5). 



The spread of the levels at each L\ scales to zero in 
the thermodynamic limit, as can be seen in Fig. [3^, and 
the spectrum has a linear shape reminiscent of a (rela- 
tivistic) CFT edge theory's energy dispersion. As such, 
we can compute the RES's "edge velocity," which is the 
overall slope of the entanglement spectrum. The cal- 
culation uses the following procedure: We partition the 
system into two hemispheres (0 = 7r/2) and focus on 
Na = N/2, so that the filling factor in A is identical to 
that of the original system. We then plot the mean value 
of for a given sector as a function of 1/N. From 
these values, we subtract the zero momentum energy ob- 
tained at the maximal value of L\. If RES mimics the 
edge mode's true energy spectrum in the thermodynam- 
ical limit, the extrapolated values of should be equal 
to 27rv{max.[L^] — L\) when TV ^ oo, where v is the edge 
mode velocity. The extrapolation, shown in Fig.[3]3, gives 
an extracted velocity of v = 2.254(3). 



A. = 1 Integer Quantum Hall 

As a first example, we consider the u = 1 IQH state. 
In the orbital basis, this state can always be written as 
a tensor product |111 . . . 11) = |11)^ |1 . . . 11)^, where 
I 111 . . . 11) means that each one of the A^^ + 1 = A or- 
bitals is occupied by exactly one electron. Thus, for any 
orbital partition, its OES always consists in a single state 
(see inset of Fig. [2| , which does not correspond to its ex- 
pected edge mode. This shortcoming is averted by the 
use of the real-space partition. As shown in Fig. [2j the 
RES exhibits a linear dispersion relation. (Note that all 
our entanglement spectra are displayed for Na = N/2, 
which is the one that contains the largest amount of in- 
formation.) As expected, the counting in each (Na^L^) 
sector is given by the counting of states for a number un- 
pinned bulk quasiholes corresponding to the same quan- 
tum numbers. This equals the number of ways to put 
Na particles in A^+i orbitals such that the sum of their 
orbital numbers equals L^. Note that the RES for the 
u = 1 state exhibits a chiral structure of levels. 



B. = 1/3 Fractional Quantum Hall 

We next consider the u = 1/3 Laughlin state (Fig.|4^). 
As expected, the counting of the RES levels is given by 
the counting of the quasiholes states of Na particles in 
A^ + 1 orbitals at momentum L\. The RES of this state 
is chiral and exhibits more pronounced non-linearities 
than the RES of the v = 1 IQH state, but should be- 
coming increasingly linear as the thermodynamical limit 
is approached, as shown in Fig. |5j Notice that similar re- 
sults have already been obtained for the OES in Ref. [371 

C. V — 2/3 Fractional Quantum Hall 

The RES is also able to exhibit the full edge mode 
structure of particle- hole conjugate states, which com- 
bine chiral and anti-chiral edge modes. In Fig. [6) we 
plot the RES of the v = 2/3 particle- hole conjugate of 
the Laughlin state. Unlike the OES (see Fig. [T]), the 
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FIG. 3: (a) Spread in RES levels vs. 1/iV for different L\ 
sectors for B = 7r/2. An even-odd particle number effect is 
observed, and the spread of the levels vanishes in the ther- 
modynamic limit, (b) Value of the zero-spread ES-levels [as 
shown in (a), the spread of the levels for each vanishes in 
the thermodynamic limit, which means that all the entangle- 
ment levels at a particular L\ collapse to the same value; this 
value is the zero-spread-ES level] extrapolated at the thermo- 
dynamic limit. Here n — max[L^] — L\. The velocity of the 
edge mode is = 2.254(3). 



structure of the RES reveals the presence of both a chi- 
ral and an anti-chiral modes. Starting from the maximal 
value of I/^, the entanglement levels first have an over- 
all anti-chiral shape with negative velocity, followed by 
an overall chiral shape with positive velocity. In finite 
size samples, the two modes cover different total ranges 
of L\ (see Fig. l6| due to the different compactification 
radii of the 11(1) CFTs in this model. At intermedi- 
ate values of L^, the two modes combine in finite size, 
and their counting matches that of two mixed chiral and 
anti-chiral branches. Even though the RES can detect 
the presence of counter-propagating modes, it is not reli- 
able in pining down the topological order of particle-hole 
conjugate FQH states. For example, given the RES of 
Fig. [6j and no other information, we could immediately 
see that the state has counter-propagating edge modes, 
but could not infer that the state is the particle-hole con- 
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FIG. 4: The RES of (a) the zy = 1/3 Laughlin state and 
(b) the = 1/3 Coulomb interaction ground state, both for 
N = 8 particles with Na = 4, on a sphere partitioned into 
two hemispheres (B — vr/2). 



jugate of the v = 1/3 Laughlin state. This is due to the 
fact the spectrum is filled with levels obtained by excit- 
ing some number of chiral and anti-chiral modes. How- 
ever, having earned the information that the state has 
counter-propagating modes, we could now plot the OES 
of the u = 2/3 state, which would give us the mirror- 
image of the u = 1/3 Laughlin OES (see Fig. [T]). The 
two pieces of information allow one to identify the state 
as the particle- hole conjugate of the v = 1/3 Laughlin 
state. 



D. ly = 1 Bosonic Moore-Read State 

We have also computed the RES for non-Abelian quan- 
tum Hall states. While the study of the = 5/2 is highly 
relevant for FQH in experimental systems, we will focus 
on the simpler analog that occurs for = 1 bosonic sys- 
tems. The fermionic u = 5/2 state involves the second 
(n = 1) Landau level, which would require to use the 
corresponding one-body wavefunctions, instead of those 
displayed in Eq. ([8|. The u = 1 bosonic case is still a 
purely lowest Landau level problem. This system is also 
relevant for rapidly rotating ultra-cold atomic gaseJ^. 
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FIG. 5: (a) The spread in ES levels vs for different L\ 
sectors, using B = 7r/2. An even-odd particle number effect 
is observed, and the spread of the levels vanishes in the ther- 
modynamic limit, (b) The zero- spread ES levels extrapolated 
at the thermodynamic limit. Here n — max[L^] — L\. The 
even-odd effect is just a consequence of N a being the integer 
part of N 12. The velocity of the edge mode \^ v — 1.41(5). 
Such value would be compatible with a rescaling of the v —\ 
edge mode velocity with a factor 1 / a/3. 




FIG. 6: The RES of the = 2/3 particle-hole conjugate of 
the Laughhn state for TV = 14, Na = 7, and 6 = 7r/2. 



For such systems, the physics is similar to that of FQH 



systems, except with bosons rather than electrons, and 
replacing the Coulomb interaction with the short range 
(two-body) ^-function interaction. There are several nu- 
merical studies that indicate this v — \ bosonic system 
is accurately described by the bosonic MR state^^ 

The RES for the bosonic MR is shown in Fig. [7^. As 
expected, its level counting matches that of the CFT 
edge mode counting. The corresponding RES for the 
delta- function interaction is depicted in Fig. [7]3. There 
is a clear entanglement gap and the low-lying branch has 
counting, shape, and values almost identical to those of 
the bosonic MR state. The structure of the spectrum 
above the gap will be analyzed in the following section. 



IV. HIGHER LEVELS IN THE 
ENTANGLEMENT SPECTRA 

Similar to the situation for the OES, the RES of the 
Coulomb interaction ground state at v = 1/3 (shown 
in Fig. [4]3) exhibits the same low-lying structure as 
the RES for the v = 1/3 Laughlin state. Additional 
branches are also clearly observed in the higher ES levels. 
These branches form what was thought to be the "non- 
universal" part of the Coulomb ES and are absent in the 
Laughlin model state ES. Similar to Ref. |42] for OES, 
we find that the RES branches of the Coulomb interac- 
tion ground-state are organized in a hierarchical struc- 
ture that mimics the excitation-energy structure of the 
model pseudopotential Hamiltonian that has its ground 
state given by the Laughlin state. These structures can 
be accurately modeled using quasihole-quasielectron ex- 
citation wavefunctions. The study of Ref. can be re- 
peated here for the RES, with similar results. 

The non-Abelian states also exhibits a series of higher 
energy branches that mimic the excitation-energy struc- 
ture of the model pseudopotential Hamiltonian that 
has its ground state given by the MR state. These 
structures can be accurately modeled using the non- 
Abelian quasihole-quasielectron wavefunctions described 
in Ref. 43 (these excitations can be obtained by tak- 
ing two layers of Laughlin v = 1/2 states, forming 
quasielectron-quasihole excitations in these layers, then 
symmetrizing over the electron coordinates between the 
layers), and the results are presented in Fig. [7| 



V. REAL SPACE ENTANGLEMENT 
SPECTRUM OF QUASIHOLE EXCITATIONS 

We analyze the change in the RES for the MR state as 
a non-Abelian quasihole is taken from the north pole to 
the south pole across the real-space partition boundary. 
The evolution of the OES under this operation was an- 
alyzed in Ref. |44| The RES transformation is similar to 
that of the OES. Partitioning the system in real space at 
O = 7r/2, for the state with two quasiholes at the north 
pole in the vacuum (/) channel (equivalently, one Abelian 
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FIG. 8: (a) The RES for the (z/ = 1 bosonic) MR state with 
N = 14 and one Abehan quasihole at the north pole. The 
wavefunction for this configuration is the Jack polynomial 
with root partition 02020202020202. (b) The RES for the 
MR state with N — 14 and one non-Abelian quasihole at the 
north pole and one non-Abelian quasihole at the south pole. 
The wavefunction for this configuration is the Jack polyno- 
mial with root partition 11111111111111. For both of these, 
the sphere was partitioned into two hemispheres (G = '7r/2) 
and the particle sector is Na = 6. 



VI. REAL SPACE ENTANGLEMENT 
SPECTRUM ON THE ANNULUS 



FIG. 7: The RES of (a) the = 1 bosonic MR state, (b) 
the ^-interaction ground state, and (c) an approximate state 
of ^-interaction ground state involving up to 2-quasihole-2- 
quasielecton excitations for N = 12 particles with Na = 6. 
The sphere was partitioned into two hemispheres (G = 7r/2). 



quasihole sits at the north pole), we see an entanglement 
spectrum that mimics that of the pure MR ground-state, 
with the counting 1,1,3... (see Fig. [8|. Taking one of 
these non-Abelian quasiholes across the sphere to the 
south pole brings us to the state with one quasihole at 
the north pole and one at the south pole. The counting 
in the RES is now 1,2,4... which is the that of the MR 
state with a a charge on the edge. The change in the 
counting occurs when the particle crosses the partition 
boundary. 



We now analyze the evolution of the RES when the ge- 
ometry of the sample is altered from the disk to the thin 
annulus. The annulus represents the geometry of the so- 
called "conformal limit" for the OES, in which one no- 
tices a complete separation of the low-lying model-state 
universal ES levels from the high-lying non-universal lev- 
els, for the Coulomb ground states of bosons at = 1/2 
and 1 and fermions at = 1/3 and 5/2. We start from 
the filled disk and then proceed to insert a large number 
of quasiholes in order to reach the thin annulus limit. For 
a state of fixed number of particles (we work with = 8) , 
this pushes the FQH liquid into an annulus of thickness 
AR = N/{27TRjii). For ^ oo, the radial spatial ex- 
tent of the annulus vanishes, the partition boundary is 
close to both edges of the annulus, and the entanglement 
spectrum is influenced by the presence of the annulus 
edges. This gives rise to the chiral and antichiral spectra 
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FIG. 9: Evolution of the RES with the changing of sample 
geometry from the disk to the thin annulus. (a) The RES for 
the Laughlin state with N — 8 for the disk geometry, with a 
partition made at radius R = y^A^^. The RES exhibits the 
usual chiral structure, (b) The RES for the Laughlin state 
on a thick annulus, obtained by inserting 20 quasiholes at the 
origin. Note that the spectrum has flattened from the case 
of the disk, (c) The RES for the Laughlin state on the thin 
annulus, obtained by inserting a infinite number of quasiholes 
at the origin. The spectrum has a flat portion, sandwiched 
between a chiral and anti-chiral region (see the blue and green 
arrows) of states corresponding to the edge modes on the two 
sides of the annulus. 
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FIG. 10: Comparison between the entanglement entropy us- 
ing the real-space partition and the orbital partition for the 
= 1/2 bosonic Laughlin state with N = 11 particles. 



VII. ENTANGLEMENT ENTROPY 



Using the real-space partition, we have computed the 
entanglement entropies for different model states in order 
to numerically verify the perimeter law and extract the 
value of the topological entanglement entropy. We com- 
pared the entanglement entropy obtained using a real- 
space partition with that obtained using an orbital par- 
tition, for the bosonic u = 1/2 Laughlin state. As can be 
seen in Fig. [Toj the difference between these can be quite 
substantial. 

We have also attempted to numerically verify the for- 
mula of Ref. 28 for the entanglement entropy of the u = 1 
IQH state in the thermodynamical limit: Sa — CyC — 
with a geometry-dependent Cy = 0.203 and 7 = 0. For 
various boundary lengths C we compute the entangle- 
ment entropy for up to = 20 particles, then calculate 
the thermodynamic value using a polynomial fit in 1/A; 
this gives the thermodynamic limit of the entanglement 
entropy for a boundary of length C The estimated val- 



ues and their errors are plotted in Fig. [TT] along with a 
linear fit, with only C > 10 values taken into account, of 
these data and the theoretical prediction for the C ^ 1 
regime. The linear fit of our data gives 7 = 0.241(5), 
which is not compatible with the expected value 7 = 0. 
However, the difference between our data and the analyt- 
ical values vanishes for large C (see inset of Fig. 11 ). This 



in Fig. [9| 



indicates that, even for this very simple and ideal case, 
we are not able to reliably extract the value of the topo- 
logical entanglement entropy without prior knowledge of 
its value (which would allow informed, but biased deci- 
sions on what data to keep) for the attainable system 
sizes in numerical studies. This casts doubt on prior nu- 
merical studies that claim to have extracted values of 7 
that match the theoretical predictions. 
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FIG. 11: Entanglement entropy as a function of the length £ 
of the boundary between the two subspaces. The red Xs are 
the (thermodynamic limit extrapolation of the) entanglement 
entropies for a given boundary length. The dotted green line 
is a linear fit of the entanglement entropies for £ > 10. The 
solid blue line is the analytical result from Ref. l28] that is 
valid in the C large limit. In the inset, we show the difference 
between the analytical and numerical values. 



VIII. CONCLUSIONS 



particle-hole conjugate states, both of which fail to show 
up in the OES. For Coulomb states at fractional filling, 
the RES exhibits an entanglement gap and a series of 
higher entanglement energy branches whose level count- 
ing is in one-to one correspondence with that of the hier- 
archical quasihole-quasielectron excitations of the model 
state. This is similar to the information given by the 
OES; however, in the latter case, a full entanglement gap 
completely separating the low entanglement energy, uni- 
versal levels from the spurious level, can sometimes exist 
by going to the conformal limit. In the case of the RES 
an entanglement gap exists only for the very few first lev- 
els of the spectrum. We have also analyzed the behavior 
of the RES as a non-Abelian quasihole is taken from the 
north pole to the south pole, passing across the partition 
boundary. The RES exhibits a change in its counting as 
the quasihole crosses the boundary, corresponding to the 
change between different topological sectors of the CFT 
(i.e. different topological charge on the boundary). Fu- 
ture research will focus on the real-space partitioning on 
the torus, on which the OES was investigated in Ref. |45] 
and the PES in Ref. |^ 

Note added: During the completion of this manuscript, 
we learned of similar unpublished works^-^*^. 



In this paper we have analyzed the entanglement spec- 
trum for a sharp, real-space partition of the space in 
which a FQH liquid exists. We have showed that the 
RES is intimately related (through diagonal transforma- 
tions involving geometry dependent normalization fac- 
tors) to the PES and that, in particular, the two en- 
tanglement spectra have identical level counting. This 
counting is bounded by the number of CFT edge mode 
levels of the particular model state analyzed. It has been 
proven that this bound is saturated in the thermody- 
namic limit for these states and also at finite sizes for 
some of these states. We expect that the bound is satu- 
rated at any finite size for all model states described by a 
CFT, and have numerically confirmed this for all sizes ex- 
amined. We have also shown that the RES offers comple- 
mentary information to the OES. In particular, the RES 
exhibits the expected (gapless) CFT edge spectra, even 
for the u = 1 filled Landau level (IQH) state, and also re- 
veals the presence of counter-propagating edge modes in 
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The total quantum dimension V = where da 

are the quantum dimensions of the anyonic charge types 
that can exist in the corresponding TQFT that describes 
the long-distance physics of the topologically ordered sys- 
tem. It is also related to the dimension of the ground- 
state subspace 1-Lg for the system on a genus g surface 
by dimTig = with 1 < Cg < V'^ for all g. For 

a Z2-graded TQFT, i.e. a topological spin theory, which 
are relevant for the description of fermionic quantum Hall 
states, the anyonic charges come in Z2 doublets (e.g. the 
vacuum and the electron charges form a doublet), so one 
should use V = as the total quantum dimension of 

the Z2 -graded theory, in order to count each Z2 anyonic 
charge doublet only once. 



